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Abstract 

We present the design, data and results from the NEXT prototype for Dou- 
ble Beta and Dark Matter (NEXT-DBDM) detector, a high-pressure gaseous 
natural xenon electroluminescent time projection chamber (TPC) that was 
built at the Lawrence Berkeley National Laboratory. It is a prototype of the 
planned NEXT-100 136 Xe neutrino-less double beta decay (Ou/3/3) experiment 
with the main objectives of demonstrating near-intrinsic energy resolution at 
energies up to 662 keV and of optimizing the NEXT-100 detector design 
and operating parameters. Energy resolutions of ~1% FWHM for 662 keV 
gamma rays were obtained at 10 and 15 atm and ~5% FWHM for 30 keV 
fluorescence xenon X-rays. These results demonstrate that 0.5% FWHM res- 
olutions for the 2,459 keV hypothetical neutrino-less double beta decay peak 
are realizable. This energy resolution is a factor 7 to 20 better than that 
of the current leading Ou/3/3 experiments using liquid xenon and thus repre- 
sents a significant advancement. We present also first results from a track 

2 



imaging system consisting of 64 silicon photo-multipliers recently installed in 
NEXT-DBDM that, along with the excellent energy resoution, demonstrates 
the key functionalities required for the NEXT- 100 0^/3/3 search. 
Keywords: Xenon, HPXe, Energy Resolution, High-Pressure, TPC, 
Electroluminescence, Double Beta Decay, Neutrinoless, 136 Xe, Fano, 
gaseous detectors for imaging 

1. Introduction 

Neutrino-less double beta decay (0vf3f3) is a postulated pQ rare process 
in which a nucleus changes by two units of charge while emitting two elec- 
trons (or positrons) without the emission of neutrinos |2j. Should this decay 
happen in nature the sum of the energies of the two electrons will be mo- 
noenergetic at exactly the Q-value of the nuclear decay (Qpp, equal to the 
mass difference between the parent and daughter nuclei). A precise energy 
measurement can therefore greatly aid in the identification of the 0^/3/3 pro- 
cess in the presence of other more common processes that produce either 
continuous energy deposition spectra or peaks at well known and well sepa- 
rated energies. The occurrence of Qv/3/3 would imply that neutrinos are their 
own anti-particle [21 H], or Majorana particles. Should neutrinos prove to be 
Majorana particles the observed prevalence of matter over anti-matter in our 
universe could be explained through the Leptogenesis mechanism (see e.g. 

mm)- 

The 136 Xe isotope is one of the candidate nuclei in which 0v(3f3, in this 
case 136 Xe — > 136 Ba + e~ + e~ , could be detectable because the single beta 
decay which would otherwise dominate the experimental count rate is en- 
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ergetically forbidden. Two currently running experiments are searching for 
this rare process in 136 Xe. EXO [7j uses 200 kg of the enriched isotope in a 
cryogenic liquid xenon TPC and KamLAND-Zen [8] uses 330 kg of the iso- 
tope dissolved in 13 tons of organic scintillator. The energy resolutions for 
these two experiments are 3-4% and 10% FWHM respectively at the 2.459 
MeV [9] Q/3/3 of the 136 Xe decay. It is well known, on the other hand, that 
xenon in gaseous phase can offer significantly better energy resolution due 
to its small Fano factor [TU] F = 0.14 (a measure of the level of fluctuations 
in the number of ionization electrons). For a xenon gas pressure of less than 
57 atm at room temperature (density of 0.55 g/cm 3 ) the intrinsic energy 
resolution is expected to be [UJ about 0.3% FWHM near the Qpp. It is thus 
clear that a xenon detector at moderately high pressure would represent a 
significant advantage for the search of the 0v/3/3 spectrum peak as long as its 
implementation can preserve a near-intrinsic energy resolution. 

NEXT-100 is an experiment [12] being constructed to search for 0^/3/3 
using 100-150 kg of 136 Xe in a 10-15 atm TPC at the Canfranc Underground 
Laboratory (LSC) under the Pyrenees mountains in Spain. In the TPC, 
conceptually developed in Ref. [32], electrons liberated through ionization by 
the passage of energetic charged particles (such as the two electrons from the 
Qv(3(3 decay) drift under the presence of a weak electric field towards a thin 
region with a high electric field. The E/P (electric field divided by pressure) 
in this high field region is such that electrons acquire enough energy to excite 
xenon atoms, but not enough to ionize them. Most of the excitation energy is 
ultimately released in the form of ultra-violet (VUV) photons of wavelengths 
near 172 nm and constitutes the electroluminescence (EL) [14J signal. For 
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each ionization electron, thousands of EL photons can be produced in the EL 
amplification region, with extremely low fluctuations in the number of EL 
photons produced. An array of photomultiplier tubes (PMTs) then detects 
a fraction of the VUV photons to render a measurement of the total energy 
released in the gas with a statistical precision near the Fano limit. 

The NEXT-100 TPC will provide, in addition to a very precise energy 
measurement, a 3-D image of the ionization tracks by means of a dense array 
of silicon photomultipliers [TJ)] (SiPMs or MPPCs) installed near the elec- 
troluminescence region. This topological information is useful to distinguish 
between events with 2 electrons emerging from a single point, such as in 
Of 0/3, from events with one or more electrons that result from single-site 
and multi-site interactions of gamma rays from natural radioactivity in the 
detector and surrounding materials. 

In this paper we present the design, data and results from the NEXT- 
DBDM (NEXT prototype for Double Beta and Dark Matter) detector, a 1 
kg natural xenon electroluminescent TPC that was built at the Lawrence 
Berkeley National Laboratory. It is a prototype of the NEXT-100 detector 
with the main objectives of demonstrating the near intrinsic energy resolution 
at energies up to 662 keV and of optimizing the NEXT-100 detector design, 
construction, and operating parameters. 

2. High Pressure Xenon Electroluminescent TPC 

In the past, various high pressure xenon detector designs have been cre- 
ated for studies directed towards the improvement of energy resolution. Stud- 
ies of gridded ionization chambers in cylindrical form [16] and with parallel 
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plate geometries [TTJ Q2J [THl E] were able to achieve energy resolutions at the 
level of 2-2.5% FWHM at 662 keV. Later studies investigated alternatives to 
shielding grids such as a hemispherical geometry [20J, using a virtual Frisch 
grid [21J, pulse correction using drift time determined by observing primary 
scintillation [22] , and variations [211 23] of the coplanar anode approach in- 
troduced in [24j. These studies produced energy resolutions of 2-6% FWHM 
at 662 keV. It has been shown in [2S] that good energy resolution (2.7% 
FWHM at 122 keV) is achievable using electroluminescent readout (see [13] 
for a more thorough historical account of electroluminescence in detectors), 
and here we extend the usage of this technique to higher energies. 

The basic building blocks of the NEXT-DBDM xenon electroluminescent 
TPC stainless steel pressure vessel, a gas system that recirculates and 

purifies the xenon at 10-15 atm, stainless steel wire meshes that establish 
high-voltage equipotential planes in the boundaries of the drift and the EL 
regions, field cages with hexagonal cross sections to establish uniform electric 
fields in those regions, an hexagonal pattern array of 19 VUV sensitive PMTs 
inside the pressure vessel and an associated readout electronics and data 
acquisition (DAQ) system. 

When ionizing radiation traverses the drift region of the TPC, xenon 
atoms are ionized or excited. Most of the excitation energy is promptly 
released as a fast scintillation pulse of 172 nm VUV photons that lasts 10-30 
ns [H]. A fraction of these photons are detected in the PMT array, forming 
the SI signal that provides the start time to for the TPC. The ionization 
(or secondary) electrons, on the other hand, drift for a maximum distance 
of 8 cm at a velocity of ~ 0.1 cm//xs towards the EL region. There, they 
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accelerate and produce copious EL VUV photons. The same PMT array 
detects a fraction of these photons, forming the S2 signal. 

In the NEXT-DBDM detector the PMT array and the EL region, which 
are both hexagonal areas with 12.4 cm between opposite sides, are 13.5 cm 
away from each other (see FigjTj) . Thus point-like isotropic light produced 
in the EL region illuminates the PMT array with little PMT-to-PMT varia- 
tion. This geometric configuration also makes the illumination pattern and 
the total light collection only very mildly dependent on the position of the 
light origin within the EL region. The diffuse reflectivity of the TPC walls 
increases this light collection uniformity further. As a result, the device pro- 
vides good energy measurements with little dependence on the position of 
the charge depositions. On the other hand, without a light sensor array near 
the EL region precise tracking information is not available and only coarse 
average position can be obtained using the PMT array light pattern. In this 



report, only the tracking analysis in section 8.1 was performed using data ac- 
quired with the tracking sensor array recently installed in the NEXT-DBDM 
TPC. In a future publication we will report on the energy resolution achieved 
in a combined analysis using information from both the PMT and tracking 
arrays. 

3. Intrinsic Energy Resolution in the HPXe TPC 

The intrinsic energy resolution in a xenon gas detector that measures 
ionization is given by: 

5E/E = 2.35^(FW i /E) (FWHM) (1) 
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Figure 1: NEXT-DBDM electroluminescent TPC configuration: An array of f9 
photomultipliers (PMTs) measures Sf primary scintillation light from the 8 cm long drift 
region and S2 light produced in the 0.5 cm electroluminescence (EL) region 13.5 cm away 
from the PMTs. Two 5 cm long buffer regions behind the EL anode mesh and between 
the PMTs and the cathode mesh grade the high voltages (up to ±17 kV) down to ground 
potential. 

where E is the energy released in the detector, W\ is the average energy 
required to liberate an electron and F is the Fano factor that quantifies the 
fluctuations in the number of liberated electrons. F and W\ are energy, drift 
field and pressure dependent due to electron-ion recombination and due to 
the energy dependence of the energy loss rate dE/dx (see for example |1 1] . [26] 
and [2 7J ) . In this study we explored the 30-662 keV energy range at pressures 
between 10 and 17 atm, and drift fields in the 0.3-2.0 kV/cm range. For the 
purpose of studying the energy resolution in the NEXT-DBDM detector we 
used 662 keV gamma rays from a 137 Cs source and 29.1-34.5 keV xenon 
X-rays |28j that follow photoelectric interactions of the gamma rays. For 
F=0.14 [29] and Wi=24.8 eV [301 H3] the intrinsic energy resolutions are 
approximately 0.53% FWHM for 662 keV and 2.5% FWHM for 30 keV. 
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4. Experimental Aspects of the Energy Resolution in the Electro- 
luminescent HPXe TPC 

4-1. Statistical 

In the NEXT-DBDM detector the number of liberated electrons is not 
measured directly. Rather, those electrons are drifted towards the EL region 
and then accelerated to produce 0(1000) VUV photons for each electron 
crossing the EL region of which O(10) are measured as photoelectrons in the 
PMT array. This gain and measurement sequence introduces fluctuations 
beyond the intrinsic Fano limit. In Ref. [31J a formalism was developed to 
calculate the energy resolution achievable in light of those additional fluctu- 
ations: 

5E/E = 2.35^((F + G)W-J E) (FWHM) (2) 

with 

G = l/r] + (1 + o-p d )/n pe (3) 

where 77 is the average number of VUV photons produced in the EL region 
per secondary electron (or optical gain), n pe is the average number of photons 
detected (as photoelectrons) per secondary electron and a^ d is the variance 
on the charge measured for single photoelectrons. r\ 

1 A thorough analysis of the statistics of energy resolution gives the formula 

( 2.^1^7^ + ^ 2 + 1 (FWHM). (4) 

\ E J V nion UEL ndct 

where n; on = E/Wi is the number of electron-ion pairs produced, tiel = "■ion?/ is the 
total number of EL photons, and n^ ct = ?yeri ion is the total number of photoelectrons 
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The electroluminescence optical gain in pure high pressure xenon is given 
approximately by [32], [13] : 



where AV is the voltage difference between the mesh planes that form 
the EL region in kilovolts, p is the pressure in atmospheres and Az is the 
distance between the electrodes in centimeters. For instance, if p=10 atm, 
AV=11.3 kV and A^=0.5 cm the EL optical gain is 77=1000. 

The value of n pe is the product of the EL optical gain 77 times the collec- 
tion efficiency (the probability of a VUV photon generated in the EL region 
reaching a PMT window) times the PMT quantum efficiency at the corre- 
sponding wavelength. For example for a 10% collection efficiency, a 15% 
PMT quantum efficiency at the 172 nm wavelength of the xenon electrolu- 
minescence and 77=1000, n pe =15 and the expected energy resolution for 662 
keV gamma rays in the absence of systematic effects is 0.66% and 3.11% for 
30 keV. 

4-2. Systematic 

Systematic effects that broaden the energy resolution can be grouped in 
two categories: position and time dependencies. Position dependencies of the 

detected by the PMTs (e is the detection efficiency of EL photons). J^p is the Conde- 
Policarpo factor which characterizes the fluctuations in EL light production similar to 
how the Fano factor F characterizes the fluctuation in ionization production. Equation [2] 
assumes Jcp = 2. In most realistic scenarios, however, the second term under the square 
root can be neglected due to the low EL photon detection probability. 
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energy response arise, for example, from a non-uniform EL light collection ef- 
ficiency, a non-uniform EL gain, secondary electron losses due to attachment 
on electronegative impurities during drift and secondary electrons hitting the 
TPC walls due to diffusion or due to drift field non-uniformities. Time depen- 
dencies of the energy response arise, for example, from time- varying voltages, 
temperature (and its subsequent effect of gas viscosity, PMT response, etc), 
gas flow and purity, PMT response and gas density. 

In the NEXT-DBDM TPC systematic dependencies of the energy re- 
sponse as a function of the position along the drift (z-axis in our chosen 
coordinate system) are small but still non-negligible due to the high xenon 
gas purity achieved. In addition, the z position of charge depositions within 
the drift region is very precisely measured. The opposite is true in the plane 
perpendicular to the drift direction (x — y): the light collection changes more 
rapidly as a function of distance from the center axis of the TPC and the 
x — y position of charge depositions is poorly measured. For this reason, to 
study the energy resolution achievable with the xenon EL TPC, calibration 
gamma rays are introduced along the center z axis through a narrow aper- 
ture collimator. Still, the actual charge depositions happen over an extended 
region in x — y (and z) due to the length of the electron tracks and due 
to the multi-site depositions from Compton scatters and from xenon X-rays 
following photoelectric absorption. 
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5. Experimental Setup 

5.1. Gas System 

The main functions of the gas system are to recirculate and purify the 
xenon at pressures up to 17 atm. A magnetically driven seal-less and oil-less 
pressure rated (95 atm) pump manufactured by Pumpworks Inc. (model 
2070) recirculates the xenon at room temperature at 5-15 standard liter per 
minute (slpm). For the total system volume of approximately 10 liter at 
10 atm pressure the pump recirculates one full volume in about 10 min- 
utes. A pressure rated (18 atm) heated getter from Johnson Matthey (model 
PureGuard) removes O2, H2O, N2 and many other impurities using a non- 
evaporable zirconium-based material. The getter operates at a constant 450 
degrees Celsius irreversibly removing the impurities through bulk diffusion. 

In Figure [2] the complete gas system is shown. Besides the recircu- 
lation pump and the heated getter the system includes a vacuum system 
with a roughing pump (make TriboVacDry), a turbomolecular pump (make 
Leybold-Heareus) and Pirani and ion vacuum gauges, a reclamation cylinder 
where xenon is stored after it is cryogenically removed from the main pres- 
sure vessel, an argon purge system, a gas sampling system with a precision 
leak valve and a residual gas analyzer (SRS model RGA100) and a room 
temperature secondary getter (SAES model MicroTorr MCP190). 

A set of pressure relief valves (with different settings for the various parts 
of the pressure rated system) and a burst disk in the vacuum system protect 
the equipment and personnel from overpressure hazards. 

The typical gas cycle during normal operation of the TPC consists of a 
vacuum step to 10~ 5 torr, followed by an argon purge and recirculation at 1 
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Figure 2: Gas system schematic (simplified): The bottom left of the diagram shows 
the main pressure vessel with its connections to the vacuum, fill, recirculation/purification, 
sampling and reclamation components. 

atm, a second vacuum step to remove the argon and then the xenon fill from 
the reclamation cylinder. After the fill, the recirculation and purification 
is started, monitored with gas flow meters (models Sierra Smart-Trak and 
Omega FM1700) and controlled with a variac that powers the recirculation 
pump to set the recirculation flow. 

5.2. Pressure Vessel and Feedthrough Ports 

The main pressure vessel is an 8.7 liter stainless steel cylindrical shell of 
20 cm diameter and 33.5 cm length with an elliptical head on one end and a 
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custom gasketed Conflat flange (main flange) on the other. Half inch VCR 
ports on the side of the vessel are used for gas fill, recirculation flow, pressure 
and temperature gauges and pressure relief valves. A 5.9 cm diameter (ID) 
port is used for vacuum pumping. On the main flange there are small ports 
with commercial high voltage feedthroughs (rated to 20 kV and 17 atm) with 
additional custom PTFE sleeves on the inside to increase the breakdown 
voltage when operating with high pressure xenon. A larger central port is 
used to connect to an octagon-shaped vessel through a long tube with an 
internal source reentrance tube with a 2 mm thick endcap. Signal and high- 
voltage coaxial cables from the (in-vessel) PMT array pass through the axial 
extension tube and connect with 32-pin feedthroughs on the side ports of the 
octagonal vessel. 

5.3. TPC 

The field configuration in the TPC is established by five stainless steel 
meshes with 88% open area at a z position of 0.5 cm (cathode buffer or PMT 
mesh), 5.5 cm (cathode or drift start mesh), 13.5 cm (field transition or EL- 
start mesh), 14.0 cm (anode or EL-end mesh) and 19.0 cm (anode buffer or 
ground mesh) from the PMT windows. Electroluminescence occurs between 
13.5 and 14.0 cm. The meshes are supported and kept tense by stainless 
steel frames made out of two parts and tensioning screws on the perimeter. 
The TPC side walls, made out of 18 individual rectangular assemblies 7.1 cm 
wide (and 5 and 8 cm length) connecting adjacent meshes (except around 
the 0.5 cm EL gap), serve the dual purpose of light cage and field cage. Each 
side wall assembly is made of a 0.6 cm thick PTFE panel and a ceramic 
support panel. The FTFE panels are bare on the side facing the active 
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Figure 3: Cutaway schematic of the TPC. 

volume and have copper stripes parallel to the mesh planes every 0.6 cm on 
the other side. The bare PTFE serves as reflector for the VUV light (a 40-50% 
Lambertian reflectivity was measured in [33]) to increase the light collection 
efficiency Adjacent copper stripes are linked with 100 MQ resistors to grade 
the potential and produce a uniform electric field. The ceramic support 
panels are connected, mechanically and electrically, to the outer perimeter 
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of the mesh support frames and to the first and last copper stripes on their 
corresponding PTFE panel. High voltage connections to establish the TPC 
fields (HHV) are made directly to the mesh frames. 

Six short PEEK rods going through holes on the mesh frames' vertices 
secure the anode and anode buffer meshes to the main vessel flange. Three 
PEEK c-clamps (on alternate sides of the hexagon) with grooves to prevent 
HV surface breakdown attach, in turn, the anode mesh to the field transition 
mesh maintaining a gas gap of 0.5 cm. Finally, six PEEK rods support the 
field transition, cathode and cathode buffer meshes as well as the 19-PMT 
array. Mechanical tolerances obtained are better than 1 mm throughout the 
TPC geometry with better tolerances on the electroluminescent gap between 
13.5 and 14.0 cm. 

In the initial implementation of the NEXT-DBDM TPC longer PEEK 
rods supported the entire assembly going through holes in the anode and 
field transition mesh frames with PEEK spacers to maintain the EL gap. 
High voltage breakdown in the form of sparks across the gap during HV 
conditioning and at random times thereafter produced conductive paths on 
the insulator surface, requiring time consuming repair. The c-clamp solu- 
tion described above made the TPC completely resilient to the unavoidable 
occasional sparking. 

Clean gas from the purification system flows through an internal tube to 
an enclosed volume behind the PMTs and reaches the active volume through 
small dedicated PEEK tubes between the PMTs to exit the TPC through a 
port on the octagonal vessel at the end of the extension tube. 

A wide range of HHV voltages were used to investigate the TPC perfor- 
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mance dependence on drift electric field and on the E/P in the EL region: 
cathode voltages from -4 to -13 kV, field transition voltages from -1 to -10 
kV and anode voltages from 1 to 13 kV. 

5.4- PMTs, Readout and Data Processing 

The PMT array is composed of nineteen 2.54 cm diameter Hamamatsu 
7378A PMTs. These PMTs, with fused silica windows, have a quantum 
efficiency of ~15% for 172 nm xenon light. The PMTs were individually 
pressure tested to 19 atm (absolute) and no mechanical or performance fail- 
ures were observed. The HV, typically about -1 kV, is individually set to 
get a ~ 10 6 gain. The bases for the PMT array are implemented in a single 
hexagonal PCB board with surface mount components and with pin sockets 
that provide both mechanical support for the PMTs as well as the necessary 
electrical connections. The base has a total resistance close to 1 Mf2, thus 
power dissipation is about 1 Watt per PMT. 1 /iF capacitors connected to 
the last dynode stages keep the gain constant during long EL light pulses. 
200 pF capacitors connected from the PMT anodes to ground stretch the 
pulse waveform such that all photoelectrons are properly sampled in the 100 
MHz digitizers. 

PMT anode signals travel through ~lm long PTFE coaxial cables to the 
32-pin feedthoughs and then to 8-channel Phillips-777 amplifiers set to a gain 
of 40. The amplified signals are then stretched and attenuated in a passive 
RCR circuit to reduce high frequency noise and to match the input range of 
the SIS3302 16-bit digitizers that sample the individual PMT signals at 100 
mega-samples per second. PMT waveform data, typically 16,384 samples or 
163.84 /is, are readout through an SIS3150 USB to VME interface to a Linux 
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server for processing, analysis and storage. The overall system noise is such 
that individual photoelectrons can be detected and the instantaneous (10 ns) 
dynamic range per PMT is l-to-200 photoelectrons. Custom-written DAQ 
software is used to control the data acquisition. 

The trigger is designed to identify S2 EL pulses which have much larger 
area than the SI pulses. The signal from a single PMT is integrated by a 
12/is peaking-time shaping amplifier and then sent to a discriminator. This 
scheme effectively provides an S2 energy threshold. The trigger is used as a 
stop signal (with a stop delay) for the digitizers such that there are, typically, 
80 /is of waveform data preceding the first S2 pulse and 80 /xs after. Since the 
maximum drift time is about 80 [is, this permits the search for the SI pulse 
in the offline analysis while ensuring that all the S2 light from one event is 
measured. 

After a block of 512 events is collected with about 100 MBytes in raw 
DAQ data format, an automated process unpacks the data, applies calibra- 
tion constants to the individual PMTs and executes the analysis code based 
on ROOT [3l] and FMWK [35] that finds SI and S2 pulses and computes 
energies, times and positions and outputs the results in 1 MBytes ROOT 
data trees. 

5.5. Controls and Monitoring 

All system controls (except the single PMT HV power supply), such as 
HHV voltage settings and current limits and recirculation flow, are done 
manually. The three HHV voltages and currents, the PMT HV power supply 
voltage, the pressure inside the TPC vessel, the temperature at two points in- 
side the TPC and the room temperature are automatically and continuously 
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recorded for monitoring and to aid in offline data analysis. The recirculation 
flow is recorded manually. 

6. Data Analysis 

In the first step of the data analysis the individual PMT waveforms are 
changed to photoelectron units using calibration constants from dedicated 
low occupancy single photoelectron runs with short LED light pulses. A 
sum waveform is then computed from a sample-by-sample addition of the 
19 PMTs' waveforms and the baseline of the sum waveform is obtained. A 
search for SI and S2 pulses that cross a threshold follows. This threshold 
is determined from the value of the baseline noise. Pulses are defined as 
SI candidates if they are less than 500 ns wide and have a maximum rise 
time (defined as the time from the start of the pulse to the location of the 
peak amplitude) of 100 ns. All other pulses are considered S2 candidates. 
Individual pulses are integrated and the largest SI candidate is assumed to 
be the event's start time and the others discarded. All S2 pulses that follow 
the chosen SI are considered part of the event while the ones preceding it 
are discarded. An event is considered valid if it has an SI pulse with at least 
one associated S2 pulse. Figures [4] and [5] show typical valid events from 662 
keV 137 Cs gamma ray interactions. 

During the automatic analysis of the data, the electron attachment life- 
time of the gas is not yet known. In order to enable the offline correction 
of these charge losses the ten first moments of the S2 charge distribution 
J q(t)t n dt with n=0,l...,9 are calculated where t is the time interval since the 
SI start time pulse. 
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Figure 4: Typical full energy 662 keV gamma ray event waveform: the sum of the 

previously calibrated waveforms of the 19 photomultiplicrs is shown. The SI pulse (shown 
in detail in the left inset) due to xenon scintillation is short and with 0(200) measured 
photons. Two S2 pulses caused by electroluminescence of xenon from ionization electrons 
follow. The first with ^270,000 measured photons is likely due to a 630 keV electron from 
the photoelectric interaction of the 662 keV gamma ray; the pulse structure reflects the 
ionization density of the track along its ^2.5 cm long projection on the drift (z) direction. 
The second (shown in detail in the right inset) with ^12,000 measured photons is likely 
due to the interaction, a few cm away, of a 30 keV xenon X-ray following the photoelectic 
process; since the 30 keV energy deposition is nearly point-like the pulse shape is gaussian 
with a <tl of 1.4 mm set by the longitudinal diffusion of the electrons during the ~8 cm 
drift. This event is from a data run taken at 10 atm with a 0.16 kV/cm drift electric field 
and an E/P of 2.1 kV/(cm atm) in the EL region. 

An average x — y position for the event is calculated from an S2-charge 
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Figure 5: A second typical full energy event: In this event from the same data run 
as in Fig. g]thc xenon X-ray interacted close to the EL region. The ionization electrons 
from it drifted for just 3 mm and underwent a correspondingly small longitudinal diffusion 
a i, of 0.3 mm. 

weighted average of the PMT positions. The event is also time split into 
equal-charge slices. For each slice an average x — y position is computed 
from the PMT light distribution. This provides x — y — z positions for well 
separated energy depositions such as from Compton scatters or from X-rays 
following photoelectric interactions. 

7. Results 

7. 1 . PMT Performance 

As shown in formulas [2] and [3] the energy resolution achievable in an elec- 
troluminescent TPC depends on the precision with which individual photons 
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can be measured. System noise and the variance from the avalanche multi- 
plication in the PMT (mostly from the first dynode stages) contribute to the 
spread of the charge measured for photoelectrons. We found, however, that 
the charge fluctuations due to PMT after-pulsing are the dominant factor 
in the photoelectron charge resolution. After-pulsing caused by the occa- 
sional ionization of residual gas molecules in the PMT vacuum volume can 
produce delayed pulses with large charges (10-20 photoelectron charges are 
not uncommon). Dedicated LED data runs with light pulses less than 100 
nsec long (first afterpulse peaks from hydrogen and helium ions appear at 
200 and 300 nsec respectively for these small PMTs) were thus taken to as- 
sess the charge-variance for each PMT. The typical value was found to be 
a{Q)/Q = 1.2 thus o^ d =1.44 with small PMT-to-PMT variations. 

7.2. Position Measurement 

The NEXT-DBDM TPC configuration with the PMT array 13.5 cm from 
the EL region does not permit detailed track reconstruction in the x — y 
plane. Still, the position reconstruction achievable allows the fiducialization 
of pulses to select events/pulses within regions of the TPC with uniform 
light collection efficiencies. The position reconstruction for isolated 50-100 
keV energy depositions shown in Fig. [6] displays the hexagonal boundary of 
the TPC. A scaling factor of ~30 and an x — y offset are needed to convert 
charge-weighted average x and y positions to true TPC (x, y) coordinates. 
The non-uniformity observed in Fig. [6] is due to the inaccuracies of the simple 
charge- weighting algorithm. A detailed Monte Carlo simulation of the exper- 
imental setup, including the EL light propagation and the optical properties 
of the detector components shows the same non-uniformities, a nearly linear 
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correspondence between true and reconstructed positions and the same ~30 
scaling factor between them. The x — y offset, on the other hand, is unique 
to the TPC data and is due to small PMT-to-PMT differences in quantum 
efficiency and gain (not fully removed by the dedicated calibration). 
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Figure 6: Position reconstruction: The charge-weighted average of the 19 PMT 
positions is used for x — y reconstruction. Events with energy depositions between 50 and 
100 keV were selected; at these energies ionization tracks extend for just a few mm and 
produce enough light to reconstruct with sufficient position resolution. The edges of the 
hexagonal area correspond to the TPC cross section. Due to the spatial uniformity of the 
PMT plane illumination a scaling of factor of ~30 (not applied here) is needed to obtain 
true x — y positions from these charge-weighted averages. These data were taken with the 
137 Cs source, at 10 atm and with E/P of 2 kV/(cm atm) in the EL region. 

Several data runs were taken with LED light pulses of various intensities 
to determine the position resolution for point depositions of charge as a 
function of the amount of light produced and detected. Figure [7] shows the 



23 



obtained position resolutions which approximately follow the expected 1/y/N 
dependency on the photon statistics. Point-like depositions of charge, such 
as those from xenon X-rays, can thus be reconstructed in x — y with better 
than 1 cm resolution if at least 10,000 EL photons are detected from them. 



E 




Figure 7: Position resolution for LED light pulses: The position resolution in the 
x — y plane as a function of the number of detected photons is shown. The resolution 
is derived from data from light pulses from a 378 nm LED located behind the PTFE 
backplane, which diffuses its light, and illuminates the PMT plane from a distance of 
about 19 cm. The inset shows the radial distributions (and fits) from which the individual 
resolution values were obtained. 

7.3. Energy Measurement 

The energy measurement is derived from the summed waveform (in pho- 
tons detected) of all S2 pulses following the event SI pulse. Figure [8] shows 
a calibrated spectrum obtained from the 137 Cs source collimated with the 
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gamma rays entering the chamber along its central axis. The full energy 
peak events in this data run have ~270,000 detected photons which, for 
W / i=24.8 eV, corresponds to about 10.1 photons detected per ionization (sec- 
ondary) electron. Using the nominal EL gain r/=846 from Eq. [5] the total 
number of photons produced is 846 • 662, 000/24.8 ~ 22.6 million. From this, 
an S2-photon detection efficiency of 1.2% is deduced. 

While the spectrum was calibrated only using the 662 keV line, the xenon 
X-rays peak appears at the correct energy (two strongest X-ray lines are at 
29.782 keV and 29.461 keV [28]) confirming the expected linearity of the 
energy measurement in the interval of interest here. 

Attachment losses of drifting secondary electrons can be assessed from the 
full energy peak position as a function of the drift time of the events. Figure 
[9] shows typical data for 10 and 15 atm with excellent electron lifetimes of 
36 and 9 milliseconds, respectively. Less than 1% of the electrons are lost 
to attachment for the maximum 8 cm drift length. Assuming that the main 
source of electron attachment is due to residual O2, the measured electron 
lifetime at 10 atm corresponds to a residual oxygen component in the 2-4 part 
per billion range |36j . The lower lifetime at 15 atm is due to the quadratic 
pressure dependence of the 3-body attachment process rate. 

The S2 photon collection efficiency is expected to vary somewhat as a 
function of the x — y position because of solid angle and chamber optics (wall 
reflectivities, etc). A sample of xenon X-ray energy depositions was used to 
measure the change in photon detection efficiency as a function of the radial 



position. Figure [10] shows this dependence derived from a run in which 
individual X-rays were reconstructed with ~0.8 cm resolution in x — y. For 
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Figure 8: Energy spectrum for 13 Cs 662 keV gamma rays: The dashed line shows 
the calibrated spectrum (using the full energy 662 keV peak alone) from a data run taken 
at 10 atm with a 0.3 kV/cm drift electric field and an E/P of 2 kV/(cm atm) in the 
EL region. The 1 mCi source is strongly collimated. Gamma rays enter the TPC along 
the drift axis (z). Below the full energy peak at 662 keV, the smaller X-ray escape peak 
at ^630 keV, the Compton edge at 477 keV and the xenon X-ray peak at ^30 keV can 
be seen. The solid line is the same spectrum with a requirement that the reconstructed 
average position of an event be less than ^1.2 cm from the TPC axis. The small feature 
near 184 keV is due to Compton backscatters. The no-source background spectrum (not 
shown) has a broad peak near 100 keV. 

the purpose of this measurement X-ray depositions can be considered point- 
like because 30 keV electron tracks at these pressures extend typically for 
less than half a millimeter |37j. A nearly flat response is seen up to 2 cm 
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Figure 9: Electron attachment losses: The average calibrated S2 charge versus the 
charge-weighted drift time is shown for full energy events. Attachment lifetimes are ob- 
tained from the exponential fits. The 10 atm data was taken with a 0.19 kV/cm field in 
the drift region and the 15 atm data with 0.59 kV/cm. 

radius followed by a linear decline in response towards the outer edge of the 
TPC where the collection is about 10% lower than in the center. 

In the absence of detailed track imaging, the radial dependence of the 
detector response to point-like depositions affects the TPC energy resolution 
for extended tracks. Electron tracks from 662 keV gamma interactions at 
the pressures of interest here have 2-4 cm spans and nearly random shapes 
due to the large coulomb multiple scattering in xenon. Figure [IT] shows the 
effect of a radial cut on the energy resolution of 662 keV depositions and 
its efficiency Since the collimated gamma rays enter the chamber along the 
main axis, events with a small average radial position have better energy 
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Figure 10: Radial dependence of energy response: Xenon X-ray energy depositions 
were selected from full-energy 662 keV events for which the first S2 pulse has energy in 
the 20-40 keV range. The reduced response at larger radial positions (also observed in a 
detailed Monte Carlo simulation of the TPC) is due to geometric and optical effects of the 
EL light collection efficiency. These data were taken at 10 atm with a 1.0 kV/cm field in 
the drift region and 2.68 kV/(cm atm) in the EL region. 



resolution primarily because they span a smaller radial range. 



Figure 12 shows the improvement in the measured energy resolution for 
662 keV depositions with increased EL light yield. As seen in Eq. [2] and [3] 
the G term dominates the resolution at low light yields (n pe small) while the 
intrinsic Fano term is the asymptotic resolution in the large light yield limit. 



In Fig. 13 the energy spectrum in the 662 keV full energy region obtained 
at 10 atm is shown. A 1.1% FWHM energy resolution was obtained for events 
reconstructed in the central 0.6 cm radius region. A small drift-time depen- 
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Figure 11: Energy resolution dependence on radial cut: Full energy 662 keV events 
are chosen and, for each, an average x — y position is computed (and scaled to true TPC 
coordinates) using all of the S2 signal. The FWHM energy resolution (data points) is 
obtained from a gaussian fit to the peak after placing a cut on the radial position (the 
abscissa). The radial cut efficiency (solid line and right ordinate axis) is also shown. These 
data were taken at 10 atm with a 0.16 kV/cm field in the drift region and 2.1 kV/(cm 
atm) in the EL region. 

dent correction for attachment losses with r = 13.9 ms was applied. The 
xenon X-ray escape peak is clearly visible ~30 keV below the main peak. 



For the spectrum taken at 15 atm, shown in Fig. 14, a 1% FWHM resolu- 
tion was obtained. At this higher pressure the xenon X-ray is less likely to 
escape from the active volume and the escape peak almost disappears. This 
resolution extrapolates to 0.52% FWHM at (3,3/3=2.459 MeV if the scaling 
follows a statistical 1 / y/~E dependence and no other systematic effect domi- 
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Figure 12: Energy resolution dependence on amount of light detected: 662 keV 

full energy peak resolutions are shown for data runs with different total S2 light yield 
(controlled primarily by the E/P in the EL region). Energy spectra were obtained by 
selecting events reconstructed in the central 1 cm radius region and with an SI signal 
with more than 100 measured photons and not more than 600. A small electron drift- 
time-dependent attachment correction is applied. The expected <^\/\N improvement is 
observed. These data were taken with drift fields of 0.5-0.6 kV/cm at 10 atm and 0.16-0.72 
kV/cm at 15 atm. 



nates. 

In order to study the EL TPC energy resolution at lower energies, full 
energy 662 keV events that had a well separated X-ray pulse reconstructed 
in the central 1.5 cm radius region were used. The energy calibration was 
done on the 662 keV full energy peak and linearity with a zero intersect was 



assumed. Figure 15 shows the energy spectrum obtained at 10 atm with a 5% 
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Figure 13: Energy resolution at 10 atm for 662 keV gamma rays: These data 
were taken at 10.1 atm with a 0.16 kV/cm field in the drift region and 2.08 kV/ (cm atm) 
in the EL region. If assumed to follow a 1/ y/~N dependence this resolution extrapolates to 
0.57% at Qf}j3=2A59 MeV. 

FWHM resolution. The spectrum was fit to the sum of four gaussians with 
relative positions and intensities fixed to the strongest xenon X-ray lines [28J . 
The absolute position of the anchor peak and the peaks' widths (all assumed 
the same) were obtained from the fit. The anchor peak is at 29.1 keV, less 
than 2% away from its nominal 29.6 keV value. 



Figure 16 summarizes our measurements and understanding of the EL 
TPC energy resolution. The lower diagonal line represents the Poisson statis- 
tical limit from the measurement of a small fraction of the photons produced 
by the EL gain while the upper diagonal line includes the degradation (mostly 
from PMT afterpulsing) due to PMT response. The circle data points show 
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Figure 14: Energy resolution at 15 atm for 662 keV gamma rays: A 1.0% 
FWHM energy resolution was obtained for events reconstructed in the central 0.75 cm 
radius region. The attachment losses correction with t = 9.0 ms was applied. A PMT 
with a clear time varying response was removed from the measurement. These data were 
taken at 15.1 atm with a 0.59 kV/cm field in the drift region and 1.87 kV/(cm atm) in 
the EL region. 



the energy resolutions obtained for dedicated LED runs with varying light 
intensities per LED pulse. The LED points follow the expected resolution 
over the two decades range studied. The two horizontal lines represent the 
xenon gas nominal intrinsic resolution for 30 and 662 keV, respectively, and 
the two curved lines are the expected EL TPC resolutions with contribu- 
tions from the intrinsic limit and the photons' measurement. Our 662 keV 
data (squares) and xenon X-ray data (triangles) taken with various EL gains 
follow the expected functional form of the resolution but are 20-30% larger 
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Figure 15: Energy resolution at 10 atm for 30 keV xenon X-rays: A 5% FWHM 
energy resolution for 30 keV was obtained These data were taken at 10.1 atm with a 1.03 
kV/cm field in the drift region and 2.68 kV/(cm atm) in the EL region. 

possibly due to the x — y response non-uniformity. Detailed track imaging 
from a dense photosensor array near the EL region, such as the one recently 



commissioned for the NEXT-DBDM prototype (see Sec. 8.1), will enable 
the application of x — y position corrections to further improve the energy 
measurement. 

In Ref. a study of the energy resolution for 662 keV gamma rays 

at pressures near condensation (30 atm and above) using a xenon ionization 
chamber found an improvement of the resolution for increased drift fields. 
Asymptotic optimum resolutions were achieved only after applying 4 kV/cm 
or larger fields. Since these large fields would be difficult to achieve in a 
detector with one meter long drift region such as the one planned for NEXT- 
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Figure 16: Energy resolution in the high-pressure xenon NEXT-DBDM elec- 
troluminescent TPC: Data points show the measured energy resolution for 662 keV 
gammas (squares), ~30 keV xenon X-rays (triangles) and LED light pulses (circles) as a 
function of the number of photons detected. The expected resolution including the in- 
trinsic Fano factor, the statistical fluctuations in the number of detected photons and the 
PMT charge measurement variance is shown for X-rays (dot dot dashed) and for 662 keV 
gammas (dot dot dot dashed). Resolutions for the 662 keV peak were obtained from 15.1 
atm data runs while X-ray resolutions we obtained from 10.1 atm runs. 



100, we investigated the resolution dependence on the drift electric field at 



10 atm. Figure 17 shows no discernible dependence in the 0.16-1.03 kV/cm 



region investigated. 
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Figure 17: Energy resolution dependence on drift field: These data were taken at 
10.1 atm with 2.1 - 2.7 kV/(cm atm) in the EL region. The energy resolution attained is 
largely independent of the drift electric field in the investigated range. 

7.4- EL TPC characterization: light yield, drift velocity and longitudinal dif- 
fusion 

While the study of the energy resolution achievable in the EL TPC was 
the main goal of this work, we pursued other measurements as well as cross 



checks against previous results by others. Figure 18 shows the linearity of 
the EL light yield as the E/P in the 0.5 cm EL region is varied from 1 to 
2 kV/(cm atm) at 15 atm and from 1 to 3 kV/(cm atm) at 10 atm. At 
these pressures xenon deviates from the ideal gas behavior at the 10% level. 
Therefore, the more appropriate variables to describe the process are density 



(N) and reduced electric field (E/N). Figure 19 shows the density normalized 



EL yield as a function of the reduced field. As expected, points from the 10 
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and 15 atm data sets follow the same linear trend with consistent slope and 
threshold. 
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Figure 18: Light yield dependence as a function of E/P in the EL region 



The electron drift velocity was obtained from the maximum drift times 
measured (which correspond to the full drift length of 8 cm) for X-ray en- 



ergy depositions. Figure 20 shows the drift velocity measured for various 
drift fields along with two Monte Carlo calculations [36J [38J using up-to-date 
xenon-electron collision cross sections. While the data points follow the trend 
in the calculations, 10% deviations are seen at the lower drift fields. Yet, the 
general agreement observed validates the pressure and field measurements as 
well as the drift field uniformity and the xenon purity. 

The longitudinal diffusion was calculated from the time spread of EL light 
of X-ray depositions (see example fits in the upper right insets of Figs. Eland 
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Figure 19: Density normalized light yield vs. EL reduced field. 



5J) and using the measured drift velocity. For X-ray pulses with long drift 
times the longitudinal diffusion is the dominant source of pulse width with 
subdominant contributions from the transit time of electrons through the EL 



gap and from the ionization track length. Figure [21] shows the dependence 
of the longitudinal diffusion on the drift field and the corresponding Monte 
Carlo calculations that mostly bracket the data. 
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Figure 20: Drift velocity dependence on drift field: These data were taken at 10.1 
atm. 
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Figure 21: Longitudinal diffusion dependence on drift field: The measured lon- 
gitudinal diffusion is presented here as the u spread along the drift direction for a 1 cm 
drift. These data were taken at 10.1 atm. 
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7.5. SI measurement 

Unlike the S2 light which is generated within a small 0.5 cm region in z, 
the SI photons are produced throughout the entire TPC volume wherever 
the primary ionization happened. As a result the SI charge (number of 
photons detected) shows large event-by-event variations. SI light collection 
efficiency is largest for events closer to the PMT array due to the larger solid 



angle. Figure |22] shows the SI signal's dependency on the average z position 
of the event. Between 0.2 and 0.45 SI photons are detected per keV of 
energy deposited. This SI yield is consistent with the known W s , the average 
energy loss required to liberate one primary scintillation photon in gaseous 
xenon [391 SOj HI] , and the photon collection and detection efficiency in our 
TPC. SI photons are produced from direct excitations of xenon atoms by the 
ionizing particle and by ion-electron recombinations. The latter component 



is, in general, field and pressure dependent. Figure 22 shows less than 10% 
differences between the SI yield of runs taken with reduced fields of 0.04 
and 0.10 kV/(cm atm) in the drift region. We conclude that, as expected, 
any recombination effects are small for electron energy depositions of interest 
here in the range of pressures and fields investigated. 



8. Further Developments and Plans 

8.1. Tracking Plane 

Since the completion of the energy resolution studies reported above we 
have finished the construction and installation of the sensor array for the 
tracking plane. It is a 1-cm-pitch 8 by 8 square pattern array of Hamamatsu 
S10362-11-025P silicon photomultipliers (SiPMs) each with 1 mm 2 active 
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Figure 22: SI signal dependence on z position: For events with full energy 662 keV 
depositions, the number of measured SI photons is shown versus the average distance of 
the ionization track from the PMT plane (obtained from the S2 charge- weighted drift time 
and the measured drift velocity). Two runs are shown with E/P in the drift region of 0.04 
and 0.10 kV/(cm atm) with very similar SI yield. 

area. The SiPMs are coated with a thin film of tetraphenyl butadiene (TPB) 
and placed in the center of the hexagonal plane 2 mm behind the electrolumi- 
nescent region of the TPC. The SiPMs produce an electronic pulse for each 
detected photon, and the TPB serves to shift the incoming 172 nm light 
from xenon to visible light detectable by the SiPMs. Details of the SiPM 
performance and coating are given in Ref. [T5]. The low power electronics 
to amplify, shape and multiplex the signals from the array was custom built 
to fit inside the pressure vessel. The digitization of the multiplexed signals 
is done outside the vessel with the same digitizers as for the PMT system. 
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Details of the operation, calibration and analysis of the SiPM array for 
the tracking plane will be reported in a future publication. Here, we present 
first results that demonstrate the system's ability to produce detailed track 



images. Figure 23 shows the relative amounts of light detected in each SiPM 
from the electroluminescence produced by an externally tagged cosmic ray 
muon that traversed the active volume of the TPC. A reconstruction based on 
a likelihood fit of the time sliced data shows the straight muon trajectory. The 
RMS error in position is estimated to be approximately 2.1 mm for individual 



reconstructed points. Figure [24] shows the z-y projection of the reconstructed 
points obtained using the time projection of the TPC. A clear 3D track is 
obtained even for this reconstruction that uses raw (not calibrated) SiPM 
signal waveforms. 



Figures 25 and 26 show the 2D and 3D projections of the track reconstruc- 
tion of a full energy 662 keV gamma ray event. The 30 keV X-ray deposition 
is well separated in z (time) and in x-y about 5 centimeters away from the 3 
cm long photoelectron track that shows the expected multiple scattering. 

Track information from this SiPM system will be used to map and correct 
the x-y energy response of the PMT plane with the goal of demonstrating a 
near-intrinsic energy resolution over a large part of the active volume of the 
detector. SiPM tracks will also be used to study the topological informa- 
tion that will enable, in the NEXT-100 detector, the identification of 0^/3/5 
candidate events while rejecting background single and multi-track events 
from gamma rays from residual radioactivity in detector and surrounding 
materials. 
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Figure 23: Cosmic ray muon track measured in the 8x8 SiPM array near the 
EL region: The 64 SiPMs covering an 8 cm by 8 cm area in the center of the TPC are 
shown as circles shaded relatively according to the total amount of light observed by each. 
The x-y coordinates at each point in the reconstructed track are given by the squares. 

8.2. Gas Additives and Nuclear Recoil Studies 

A research program is underway to study gas additives to optimize the 
detector performance. We are investigating the properties of high pressure 
xenon mixtures with CF4, CH 4 and TMA (trimethylamine). These stud- 
ies focus on diffusion (for track imaging), electroluminescence and primary 
scintillation yield, and energy resolution and will be reported elsewhere. 

We are also investigating the high pressure EL TPC response to low en- 
ergy nuclear recoils. O(10 keV) nuclear recoils are expected from WIMP 
(Weakly Interacting Massive Particle) dark matter elastic interactions with 
the xenon target. These studies, using MeV neutron sources, aim to charac- 
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Figure 24: Projection in zy of the cosmic muon track: The z coordinate is calculated 
from the time delay between the SiPM signals and the event Sf with the known drift 
velocity. 




Figure 25: Gamma ray event with photoelectron track and X-ray. 

terize the S2/S1 discrimination between electron and nuclear recoils and to 
optimize the light collection from the feeble nuclear recoil SI signals. 
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Figure 26: 3D reconstruction of gamma ray event with photoelectron track 
and X-ray. 

9. Conclusions 

We presented the design, data and results from the NEXT-DBDM high- 
pressure gaseous natural xenon electroluminescent TPC that was built at the 
Lawrence Berkeley National Laboratory. Energy resolutions of 1% FWHM 
for 662 keV gamma rays were obtained at 10 and 15 atm and 5% FWHM 
for 30 keV fluorescence xenon X-rays. The main contributions to the en- 
ergy resolution were studied and are well understood and thus justify the 
extrapolation to a 0.5% FWHM resolution for the 2,459 keV hypothetical 
neutrino-less double beta decay peak. This energy resolution is a factor 7 
to 20 better than that of the current leading 0vf3f3 experiments using liquid 
xenon. We presented also first results from a track imaging system consisting 
of 64 silicon photo- multipliers recently installed in NEXT-DBDM that, along 
with the excellent energy resoution, demonstrates the key functionalities re- 
quired for the NEXT-100 Ou/3/3 search. The construction of the NEXT-100 
detector with 90 kg of 136 Xe detector has began and first data is expected 
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for 2015. 
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